Colloidal Gelation with Variable Attraction Energy 
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We present an approximation scheme to the master kinetic equations for aggregation and gelation 
with thermal breakup in colloidal systems with variable attraction energy. With the cluster fractal 
dimension df as the only phenomenological parameter, rich physical behavior is predicted. The 
viscosity, the gelation time and the cluster size are predicted in closed form analytically as a function 
of time, initial volume fraction and attraction energy by combining the reversible clustering kinetics 
with an approximate hydrodynamic model. The fractal dimension df modulates the time evolution 
of cluster size, lag time and gelation time and of the viscosity. The gelation transition is strongly 
nonequilibrium and time-dependent in the unstable region of the state diagram of colloids where the 
association rate is larger than the dissociation rate. Only upon approaching conditions where the 
initial association and the dissociation rates are comparable for all species (which is a condition for 
the detailed balance to be satisfied) aggregation can occur with df = 3. In this limit, homogeneous 
nucleation followed by Lifshitz-Slyozov coarsening is recovered. In this limited region of the state 
diagram the macroscopic gelation process is likely to be driven by large spontaneous fluctuations 
associated with spinodal decomposition. 

PACS numbers: 



I. INTRODUCTION 

Colloidal suspensions gel if there is an attractive inter- 
actions of sufficient strength between the particles. This 
gelation transition has been the focus of intense research 
during the last decade since it plays an important role 
in many practical applications as e.g. processing of poly- 
mers or food technology. In spite of intensive effortsMSl 
in the past aimed at clarifying the nature of the gelation 
transition, the basic mechanism by which a fluid colloidal 
suspension turns into solid remains unclear. Many nu- 
merical studies have been proposed over the last decades 
which have brought a wealth of phenomenological infor- 
mation about the connection between microscopic attrac- 
tion and the gelation procesa^Hia However, analytical 
models are lacking, and therefore it is difficult to eluci- 
date the basic mechanisms and to extract scaling laws in 
analytical form. 

Some time ago, the idea has been proposed^ that the 
gelation transition may be interpreted as a "renormal- 
ized" glass transition where the growing colloidal clusters 
occupy an increasingly larger volume fraction up to the 
point at which their motions become governed by glassy 
correlation, the clusters become caged by their neighbors 
and the system becomes solid by interconnection or ran- 
dom packing of clusters. This scenario is different from 
what one observes in chemical gels where the bonding is 
permanent (in contrast with colloidal bonds that can be 
broken up by thermal energy) and percolation provides 
an excellent description of chemical gelation^. With col- 
loidal gels, however, simulations^ have established that 
the dynamics is strikingly different from that of chemi- 
cal gels and colloidal gelation cannot be understood with 



percolation concepts alone. The concept of colloidal gela- 
tion as a cluster-jamming transition has clearly brought 
progress in the modeling of the static structure-elasticity 
relation of dense colloidal gela^. However, it has not 
been implemented in an analytical model of the gela- 
tion transition that can be tested in comparison with 
experiments. The main problem resides in the diffi- 
culty of bridging the macroscopic mechanical response 
(the viscosity) with the mesoscopic level of the clusters 
and ultimately with the underlying microscopic associa- 
tion/dissociation kinetics of individual colloidal particles. 

Here we present an analytical model of colloidal gela- 
tion with variable attraction energy. The model pro- 
vides a framework which connects the level of the pair- 
attraction energy V with the mesoscopic level of the 
clusters and finally with the macroscopic mechanical re- 
sponse. Analytical laws can be extracted for the viscosity, 
the gelation time and the cluster size. These laws pro- 
vide a theoretical explanation to several observations in 
the past for which no theoretical description is available. 

This investigation was prompted by our finding that 
well-defined attractive interaction can be induced in sus- 
pensions of thermosensitive microgels by raising the tem- 
perature above the volume transition of these systemspS 
Figure [T] shows these particles in a schematic fash- 
ion: A dense network of the thermosensitive poly(N- 
isopropylacrylamide) (PNIPAM) is grafted to solid 
polystyrene core. Immersed in cold water the shell of 
these particles will swell. Above the volume transition at 
32-33°C most of the water will be expelled from the net- 
work because water becomes a poor solvent under these 
conditions. As a consequence, the shrunken shells will 
become mutually attractive and the strength of the en- 
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the master kinetic equation is used to obtain analytical 
expressions of the time-dependent cluster size distribu- 
tion and of the time-dependent volume fraction occupied 
by the clusters in the system, (v) Since the attraction is 
short-range and the hydrodynamics is screened from the 
clusters interior, clusters are assumed to behave hydro- 
dynamically like hard-spheres and the time evolution vis- 
cosity of the system is calculated using the time-evolution 
of the cluster population as input to the hydrodynamic 
description, (vi) The gelation time at which the gelation 
(fluid-solid) transition occurs is calculated analytically 
as the time at which the clusters connect into a random 
close packing and the low-shear viscosity diverges. 



FIG. 1: (color online). Schematic representation of the ther- 
moresponsive nanoparticles used in the gelation experiments. 
The dependence of the attraction on temperature (Eq. (31)) 
is also shown in comparison with the experimental data of 
RefpH V°° = 12fc_gT is the highest attraction energy reach- 
able with this system. 



suing attractive interaction can be adjusted precisely.^ 
The advantages of this system for the study of gelation 
arc at hand: Well-defined short-range attraction can be 
induced without adding an additional component by just 
raising the temperature. Moreover, the fluid state can 
be recovered by simply lowering the temperature below 
32°C. The sharp transition observed experimentally can 
be traced back to the predicted strongly nonlinear de- 
pendence of the viscosity on V. 

This article is structured as follows. First we present 
the analytical theory of colloidal aggregation and gela- 
tion with variable attraction energy. Then we discuss 
the predictions in terms of scaling laws and a state dia- 
gram extracted from the theory. Then we describe the 
gelation experiments with thermosensitive colloidal par- 
ticles. Finally, we present the comparison between theory 
and experiments. 



II. KINETICS OF AGGREGATION AND 

GELATION WITH VARIABLE ATTRACTION 

ENERGY 

A. Assumptions and steps in the derivation 

The model is based upon the following steps and as- 
sumptions, (i) Any two colloidal particles interact via a 
rectangular- well attractive potential of width S and depth 
V. (ii) The clustering process is described by a master 
kinetic equation with an effective association rate which 
accounts for bond dissociation, (hi) The association and 
dissociation rates between two particles are evaluated 
from steady-state solutions to the Smoluchowski diffusion 
equations for the two-body dynamics in a mutual attrac- 
tion potential, (iv) The so obtained analytical solution of 



B. Derivation 

1. Master kinetic equation for aggregation with 
reversible bonds 

Let us consider the master kinetic equation which gov- 
erns the time evolution of the concentrations of clusters 
of any size present in the system as a result of the micro- 
scopic two-body association and dissociation processes: 



dNk = I 
dt 2 
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where iVj is the number concentration of aggregates with 
i particles in each of them. Kf- is the rate of association 
between two aggregates, one with mass i and the other 
with mass j, while K~j is the rate of dissociation of a j + i 
aggregate into two aggregates i and j. The first term ex- 
presses the "birth" of clusters with mass k, the second 
expresses the "death" of clusters with mass k due to ag- 
gregation with another aggregate. The last two terms 
express the "death" and "birth" of k- aggregates due to 
aggregate breakup, respectively. Instead of considering 
the two dissociation terms in the master equation explic- 
itly, we can account for dissociation in an effective way by 
replacing the association constant with an effective size- 
independent rate constant K e s and dropping the breakup 
terms in the master equation. If association is controlled 
by diffusion, as we are going to see in the next section, 
the rate of association is in good approximation indepen- 
dent of the sizes of the two colliding clusters. Further, 
we also assume that dissociation is also independent of 
the clusters. These simplifications are indeed justified 
if we assume that two clusters aggregate by forming a 
bond between two particles protruding on the respective 
surfaces, such that the association/dissociation kinetics 
between any two clusters can be effectively described by 
means of K ff. The new master equation under these 



simplifications reads as: 
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Upon discrete-Laplace transforming this equation 26 , the 
analytical solution for the time evolution of the cluster 
mass distribution (CMD) reads: 



N k = 
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N denotes the number per unit volume of monomer par- 
ticles in the colloidal sol at t = 0. 9 is the characteristic 
aggregation time or lag time and is equal to: 



(4) 
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During this lag time, aggregation is slow because of bond 
breakage, and the formation of large clusters is unfa- 
vorable because N k ~ (t/8) k ~ 1 with t/9 <C 1. Hence, 
for t/9 < 1 we can safely assume that the system is 
mainly composed of monomers and of doublets or col- 
loidal dimers. If the bond dissociation process is stochas- 
tic, as we can anticipate, then after a time longer than 
9 the formation of stable bonds is possible. This has to 
be interpreted in a stochastic sense as the probability of 
large fluctuations around the average dissociation rate 
increases with time thus making possible the stochastic 
formation of long-lived bonds over a long time. 

In the next section we derive an analytical expression 
for 9 exploiting the fact that for t < 9 only monomers 
and dimers are present in the system. 



2. Effective association rate accounting for 
dissociation 

We start by considering the kinetics of reversible asso- 
ciation between two colloidal particles to form a dimer 



monomer + monomer ^ dimer. 



(5) 



The association rate be denoted by k + and the dissocia- 
tion rate by fc_. If we denote with n(t) the concentration 
of monomers at time t and with TV the total concentra- 
tion of monomers at t = 0, the evolution of n is governed 
by: 



^ = -k + n{tf - \k-n{t) 



1 



k_N 



(6) 



where we made use of the conservation condition: n 2 (t) = 
(N — n(t))/2, with n2 the concentration of dimers. With 
the initial condition n(0) = N, Eq.(6) has the following 
solution: 



n(t) 



fc_ 



2k + 2fc+ 



tanh(V^/2) + B/VA 
1 + (B/VA) tanh(V 'At /2) 



(7) 



with A = A:_(A:_ + 4k + N) and B = fc_ + 2k + N. One 
should note that k_ has dimensions of an inverse time, 
while k + has dimensions of [volume/time] because it is 
the rate constant of a bimolecular second-order reaction, 
whereas fc_ is the rate constant of a unimolccular or first- 
order reaction. 

With only monomers and dimers as for t < 9, the tem- 
poral evolution of the average hydrodynamic radius of 
the system a(t) as measured by dynamic light scattering 
(DLS) is given by^ 
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where I\ is the intensity of light scattered by a population 
of pure monomers of radius aj., and 7 2 is the intensity of 
light scattered by a population of pure dimers having a 
hydrodynamic radius a 2 . a(t) can be calculated by sub- 
stituting Eq. (7) together with the conservation relation 
n 2 {t) = (N- n(t))/2 into Eq.(8). To first order in t, the 
resulting expression reads as: 



a(t) oc 



8(iVfc+) 3 (l 2 // 1 )ai(l-ai/a 2 ) 



[2fc_ 



4/c+A- 



fc_(J 2 /Ji)(oi/o 2 )] : 



-t. 



(9) 



Upon taking the derivative and rearranging terms we ob- 
tain the standard form 



1 da(t) _ h 
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The truncation to first-order in time implies that we 
are neglecting the equilibrium plateau that ultimately is 
reached according to the law of mass action. Rigorously, 
this approximation is valid for k+N > 16fc_ as discussed 
in the Appendix A. While keeping this in mind, it is in- 
structive to consider its predictions also outside the rig- 
orous regime of validity. By comparing the previous two 
expressions, we are now able to obtain the effective asso- 
ciation rate accounting for bond-dissociation: 
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where a — (J 2 //i)(ai/a 2 ). Since 2i/I 2 = 1 + 
sin(2a 1 g)/2a 1 g, at the particle length scale one has 
I2/I1 — 1, while at the same time a\ja% — 1.38 for 
spheres. Hence, in good approximation and to make our 
formulae more transparent, in our model we take a = 1, 
which is not going to change our results neither quali- 
tatively nor quantitatively. With this result we can also 
identify the lag time: 
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This framework allows us to account for dissociation in 
the kinetics of aggregation between colloidal particles, in 
an effective way. It contains indeed both the association 
rate k + and the dissociation rate k-. 



These rates can be estimated by solving the station- 
ary equation of diffusion for the two particles in the 
frame of one of the two taken to be the origin. In the 
case of association, upon assuming stick-upon-contact as 
for short-range attraction one has: V 2 n = with the 
boundary conditions n = at r = 2a\ and n — const 
at r — > co. At steady-state, the solution for the rate 
of collision per unit volume or flux J follows upon in- 
tegration as: J = 47r(2oi)(2£))n 2 , which, upon using 
the Stokes-Einstein relation, leads to the Smoluchowski 
rate k + — (8/3)/cbT//j. The association rate is therefore 
independent of the size of the two particles or clusters 
that aggregate. This fact implies that larger particles (or 
clusters) aggregate at the same rate as smaller particles 
because the lowering of the diffusivity brought about by 
the larger size is exactly compensated by the increase in 
the collisional cross-section. When the attraction range S 
cannot be neglected, one has to solve the diffusion equa- 
tion for two particles in the field of force of a rectangular 
well of depth V and width S. The result ia^S. 
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where D = k B T /Sit^ai is the mutual diffusion coefficient 
of the particles. For short-range potentials 5 <§C 2a\ one 
recovers the Smoluchowski rate which we are going to use 
throughout this work. We should also mention that hy- 
drodynamic interactions and elastic deformation effects 
of polymer-functionalized surfaces might play a role as 
well in the very short ranged limit. Since we cannot ac- 
curately model the latter effect we choose here to use 
the classic Smoluchowski rate theory where the slowing 
down of the rate brought about by hydrodynamics near 
the surface cancels, approximately, with the speeding up 
brought about by the finiteness of the attraction range. 
For a discussion of this effect see RefPn 

In a similar fashion, the dissociation rate fc_ can be es- 
timated by solving the steady diffusion equation for the 
two bonded particles in the field of the attraction poten- 
tial. At steady-state one calculates the (Kramers) rate 
of escape of one particle from the attraction well which 
drives the dissociation event. What changes with respect 
to the association process is obviously the boundary con- 
dition, since now it is assumed that the two particles start 
in a quasi-equilibrium steady state in the attraction well. 
In good approximation, the result for the average disso- 
ciation rate is: 



fc_ = (D/S 2 )e 



-V/k B T 



(14) 



In Ref.-2SJ a theoretical justification for this formula and 
its derivation can be found. Here we chose this simpli- 
fied form and neglect non-essential prefactors in order to 
be consistent with our previous characterization of the 
colloidal systems under study-23. Eq.(14) is not accurate, 
however, as soon as one deals with shallow attraction 
energies V ~ k B T. Even in the limit V = this for- 
mula predicts a finite dissociation rate fc_ > although 
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FIG. 2: (color online). Comparison between the Kramers 
dissociation rate (continuous line) given by Eq.14 and the 
interpolation formula Eq. 15 which interpolates between the 
Kramers formula and the V — > limit. 



it is clear that no bond can be present to start with be- 
cause the particles are hard-spheres and dissociation is 
therefore instantaneous as the spheres collide. This ar- 
tifact is due to the approximation under which Eq.(14) 
is derived, whereby Kramers assumed in its derivation 
that the attraction well must be significant in order for 
the two particles to be in a quasi-equilibrium steady- 
state in the welPn To overcome this problem rigor- 
ously, one should solve the full time-dependent diffusion 
equation which however would undermine the analytic- 
ity of our approach. Hence here we propose the following 
semi-empirical formula which interpolates between the 
Kramers formula for V 3> k B T and the limit fc_ = at 
V = 0: 



fc_ = (D/S 2 )e- v/kBT + A/{V/k B Tf. (15) 

In our calculations below, we are going to use A = 10 6 s _1 
and f3 = 20, With this choice, the dissociation rate is 
equal to the Kramers formula for all attractions down to 
V/k B T > 2 and below it rapidly increases and diverges 
at V — 0. This interpolation formula is plotted in Fig. 2 
together with the Kramers formula for comparison. 

With these identifications, the expression of the lag 
time, with the explicit dependence on T, is given by: 

= {{fk B T/a)N + (D/6 2 )e- v / k ° T + A(V/k B T)-^} 2 
8[(§fc B T/ M )iV]3 

(16) 
The lag time is thus a function of the competition be- 
tween microscopic association and dissociation kinetics. 
In the limit controlled by association fc_ -c ik+N, the 
lag time is set by the time of diffusive transport as for 
irreversible diffusion-limited aggregation: 9 = 2/(k+N). 
In the opposite limit where dissociation is controlling, 
fc_ > 4k + N, the lag time goes as = k 2 _/(2k + N) 3 — 



D 2 e 2V / k sT _ finally^ w hen the condition k + N = k-, is can be readily calculated: 



exactly satisfied, which fixes V, such that the initial in- 
dividual frequencies of the forward process (association) 
and of the backward one (dissociation) are equal, the lag 
time scales as: 



6 ~ 25/(8/c + 7V) - n/(k B TN), 



(17) 



and it is inversely proportional to T through the inverse 
of the Smoluchowski aggregation rate. This is a physi- 
cally meaningful outcome because in this regime an in- 
crease of T causes the speeding up of the diffusive trans- 
port which reduces the lag time. The physical meaning of 
this result is that, in the regime of equilibrium aggrega- 
tion, the kinetics is controlled by the activated stochastic 
jump of the particles out of the attractive well which is 
the kinetically limiting process. Upon reducing the at- 
traction, the lag time increases because the formation 
of bonds requires stochastically a longer time. Vicev- 
ersa, upon increasing the attraction the lag time gets 
reduced because it becomes stochastically more likely to 
form long lived bonds on shorter times. Although it is 
tempting to identify the condition k + N = fc_ with the 
condition of detailed balance straightaway, this is not a 
rigorous identification. As discussed in the Appendix A 
where we refer to Tolman's definition^, this condition 
rather corresponds to the microscopic reversibility in the 
early time limit. 

Finally, it should be noted that, once that t > 9, 
the aggregation kinetics enters an extremely fast regime 
which is independent of the lag time and hence is the 
same for all (finite) attractions. During this fast regime 
the kinetics is relatively insensitive to the microscopic 
details of cluster aggregation and it is very difficult to 
detect the effect of the microscopic details in the macro- 
scopic properties which evolve very rapidly towards the 
solid state. 



3. Analytical solution for the clustering kinetics with 
reversible bonds 



We now have a connection between the mesoscopic 
time-evolution of the clusters and the microscopic in- 
teractions between colloidal particles which can also be 
allowed to vary with time. The CMD can be used to 
derive the time-evolution of quantities such as the aver- 
age cluster size R and the effective volume fraction in 
the system cf> defined as the fraction of volume occupied 
by the clusters at time t. In general, the radius of the 
sphere enclosing the cluster is given by: R k — a\k 1 ' d t 
where df is the cluster fractal dimension or mass-scaling 
exponent: k — q(R k /a) if , where q is a prefactor of order 
unity. Using Eq. (3) for N k , the average cluster radius 



R(t)= k =' 
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Here, Li n (x) denotes the polylogarithrrPH of order n of 
the variable x. Closed form expressions for the poly log- 
arithm of negative order are known only for the special 
case Li_i(x) = x/(l — x) 2 which corresponds to df = 1. 
Using this form we have that (0/t)Li_i(^£g) = 1 + (t/0), 
which is a significant simplification. Extrapolating this 
result to higher values of df we obtain a good closed-form 
approximat ion : 



R (t) = a It ) Li-i/df 
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(19) 



One can check that this approximation is reasonable 
quantitatively and it is always qualitatively correct in the 
whole range of df and 9. This formula implies that in the 
aggregation process there is a lag time of order 9, during 
which aggregation events are stochastically rare because 
of the low rate of successful collisions leading to bond- 
formation, where successful collisions are those which do 
not result in immediate thermally-activated bond rup- 
ture. If 9 is a small number, meaning that the lag time 
is short as for irreversible aggregation, then the kinetics 
transitions after short time to the R ~ (t/9) 1 ^ df growth 
law that has been reported experimentally for irreversible 
colloidal aggregation in the pastP^. In the limit df = 3 
of equilibrium aggregation where microscopic reversibil- 
ity is satisfied (see section IV. a and the Appendices for 
the connection between thermodynamic equilibrium and 
df = 3), by combining the above expression with Eq.(16), 
this treatment gives: 



i?(i)~a[l+<5 2/3 0D*) 1/3 ]. 



(20) 



which in the asymptotic limit correctly recovers the well 
known Lifshitz-Slyozov^l scaling R ~ i 1 / 3 for the growth 
rate in the coalescence (coarsening) regime of phase sepa- 
ration following nucleation under equilibrium conditions. 
The link between nucleation and phase separation is dis- 
cussed more in detail in section IV. a. Hence Eq.(19) is 
important because it covers all limits of colloidal aggrega- 
tion kinetics, from irreversible aggregation to nucleation 
at equilibrium, and provides theoretical justification to 
many experimental observations in the past. 

Similarly, the effective cluster volume fraction is given 
by: 
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This definition of the effective cluster volume f ractio n is 
the most used in the colloidal gelation literatur e! 35 * 36 ! An 



alternative would be the mass-weighted average volume 
fraction which is obtained by introducing a factor k in the 
sum in the numerator and is more used in the context of 
polymer gelation where the molecular weight is the key 
parameter in the gelation process. The mass-averaging 
is however less used with colloidal gelation as here the 
key parameter which controls the gelation process is the 
cluster size (and, of course, its cube which is the volume) 
rather than the mass. Using V = vN/4>, where <j> is the 
volume fraction of the colloidal gas and v the volume of 
a single particle, and Eq.Q for the CMD, we obtain: 



4> = </>(0/t)U_ 3/df 



t 



t + < 



(22) 



Eq.([22j) gives the effective cluster volume fraction as a 
function of the time-dependent interaction (accounting 
also for dissociation) embedded in the characteristic ag- 
gregation time r. The polylogarithm of order —3/df 
with 1 < df < 3 can be very accurately approximated as 
U_ 3/df (x) » x(x + l)- 1 +( 3 / d /)/(l - a0-i-(3/d/). Then 
the volume fraction occupied by clusters after some ma- 
nipulation becomes: 

4> = 4>[l + (t/9)][l + 2(t/9)]^- d f^ d f. (23) 



4- Linking the clustering kinetics with the 
macroscopic viscosity 

Consistent with our main approximation of treating 
the clusters as renormalized spheres occupying an effec- 
tive rescaled volume fraction <j>, we now describe the ef- 
fective viscosity of the system as a function of cf). This 
treatment applies to fractal clusters as well in the regime 
df ~ 2 where the screening of the hydrodynamic interac- 
tions from the interior makes them behave like effective 
spheres or spheroids, which is a known facfr^. The viscos- 
ity of the system can be estimated by treating the clus- 
ters as effective hard spheres since their hydrodynamic 
behavior is very close to that of hard spheres even for 
fractal clusters and the short-range attraction has little 
effect on the hydrodynamic viscous dissipation. Under 
these assumptions, it is possible to obtain analytical ex- 
pressions in closed for for the viscosity over the entire <fi 
range up to the cluster close packing where the system 
arrests. This approach has the advantage of providing 
analytical scaling laws for the gelation time as a function 
of the controlling parameters. 

The differential effective medium theory allows us to 
calculate the viscosity of a dense suspension starting from 
Einstein's method for calculating the viscosity of a dilute 
suspension of spheres. Because of the assumption of di- 
lute and non-interacting particles, one first obtains a lin- 
ear dependence, i.e. the Einstein formula r\ = ju(l + |<£) 
which accounts for the hydrodynamic dissipation of a 
single cluster*"". Upon introducing a small increment of 
particles in the system and accounting for their mutual 



correlations, it is possible to account for the many-body 
hydrodynamic interactions as well as for excluded vol- 
ume effects in an effective way, leading to the following 
expression^ 



-5/2 
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(24) 



This equation is a key result of this work. Here, <fi c ~ 0.64 
is the random close packing fraction of spheres at which 
the viscosity becomes infinite. The latter value is the ob- 
ject of many detailed studies aiming at its precise defini- 
tion. In particular, this value can vary depending on the 
size polydispersity 4 ^ and on the particle interactions 42 , 
in the range 0.56 -j- 0.67. For our scope these differ- 
ences are irrelevant and we have checked that they do 
not minimally affect the qualitative predictions of our 
model. Hence, consistent with our hard-sphere approx- 
imation in the viscosity calculation, we take </> c = 0.64. 
Also, one should note that upon approaching the regime 
cj) > 0.5, the system undergoes a glassy dynamical arrest 
where the clusters become caged by their neighboring 
clusters. Within this regime the viscosity still increases, 
as a power-law of 4> according to Mode-Coupling theo- 
ries 4 ^ and with an exponential dependence according to 



Adam-Gibbs theories 4 ^, before diverging at the random 
close packing. These scenarios could be implemented in 
our frameworks, in future work, to provide a more de- 
tailed description in the glassy regime. Here we are in- 
terested in the overall kinetics of the process and defer 
this type of detailed analysis of the glassy regime to fu- 
ture work. 

Upon replacing this expression in the above expression 
for the viscosity, we obtain a closed-form expression for 
the viscosity as a function of the microscopic interaction, 
of the volume fraction, and of time: 



r/ = fi 



with 



<p[l + (t/g)][l + 2(t/g)]( 8 - d /)/ d / 
1 - c 4>[l + (t/0)][l + 2(t/0)\&- d t)l d t 



c = (l - 4> c )/4>c 



-5/2 



(25) 



(26) 



Predictions of Eq.(25) for the viscosity as a function 
of time are shown in Fig. 3 for different values of the at- 
traction energy. It is seen how the viscosity diverges at a 
well defined gelation time t g , which is the time at which 
the effective volume fraction occupied by clusters reaches 
the random close packing fraction <p c . For shallow attrac- 
tions, it is evident that the viscosity remains constant 
and not too much above the value of the initial colloidal 
gas for a long time before it starts to increase. This 
observation has very important implications. In most 
experimental studies in the past the separation between 
gelling and non-gelling systems has been defined some- 
what arbitrarily. Indeed it is customary in the experi- 
mental practice to establish that a system does not gel 
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FIG. 3: (color online). Low shear viscosity calculated as a 
function of time with varying attraction energy V according 
to Eq. 25 and using Eq. 23 for the cluster packing fraction for 
the case df = 2. The existence of a lag time for all attractions, 
after which the aggregation is very fast for all attractions, is 
evident, ai = 85nm and <j> = 0.126. From right to left: 
V = lksT, V = 1.2/cbT, V = 2k B T, V = 3k B T, V = Wk B T. 



under certain conditions if it remains in a fluid state and 
does not aggregate significantly over a chosen period of 
time. Clearly, in this way the choice of the time span 
might be such that the observation time is shorter than 
the lag time, i.e. t t, s < and states that would gel after 
a time ~ 9 might be improperly classified as non-gelling. 
The expressions reported here can therefore be of help to 
experimentalists in establishing a more rigorous criterion 
for drawing experimental state diagrams of colloids. 



III. GELATION TIME 

We define the gelation time as the time at which 
r\ = oo, where r\ is the zero frequency viscosity. Of course, 
there are many criteria that have been proposed in or- 
der to identify the gelation point. One that is widely 
used is the Winter-Chambon rheological criterion based 
on the elastic moduli^, which requires a thorough char- 
acterization of the frequency response. Since here we are 
interested in the low-frequency behavior and will present 
experiments done in this limit, the 77 = 00 is a strin- 
gent criterion to identify the liquid to solid transition. 
In the present context, since we are focusing on the low- 
frequency behavior and we can deal with the viscosity in 
closed form, we choose the r\ — 00 criterion which im- 
plies that at the gelation point the system is connected 
and does not flow at least over a long time scale. 

From Eq.(25) it is possible to determine the gelation 
time t g as a function of the other microscopic parameters, 
such as the attraction, the aggregation rate and the vol- 
ume fraction. By setting the content of the main bracket 
in Eq.(25) equal to zero, in theory one could solve for 
the gelation time t g for an arbitrary df. In practice, the 



equation cannot be solved analytically for an arbitrary 
df. By studying the solution for a few special cases such 
as df — 3 and df — 2, and dropping high-order terms in 
<f>, we obtain the following approximate formula: 



e 



2[{l + c)(/)/2] d f/ 3 



(27) 



Hence the T dependence of t g is the same as for given 
in Eq. (16). As one could expect, also from the con- 
sideration of Fig. 3, the gelation time is proportional to 
the lag time 0. Using Eq. (12), we can relate the gela- 
tion time to the microscopic association and dissociation 
processes: 



t 9 = 



(4fc + (0A>)+fc-) 2 

8(fc+(0/tt))3 

2[(l + c)</>/2] d // 3, 



(28) 



Here we wrote N — cf)/v to make the full dependence 
on (j> explicit. Let us discuss the various limits of this 
formula as a function of the attraction energy first. When 
V/kgT 3> 1 the formed bonds are frozen in for a long 
time and fc_ ^C ik+(<p/v). In this limit we get a gelation 
which depends only on the diffusive transport and on 4>: 



1 



tg 16*+(0M(1 



c)(j)/2] d f/ 3 



-l-(d//3) 



In the opposite limit of shallow attraction, fc_ 
4k + ((j)/v), we get: 



(29) 
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16(fc + ((/>/v))3[(l + c)^/2] d // 3 



k 2 _ 



-3-(d f /3) 



(30) 

Finally, when the early-time microscopic reversibility 
condition is satisfied, we have that the gelation time fol- 
lows an Arrhenius law as a function of the attraction: 
t g ~ e ~v[k B 7 \ This law implies that the gelation time, 
in this regime, increases upon increasing T because the 
bonds become more short-lived at higher T which slows 
down the aggregation process. Therefore, it is clear that 
the temperature affects the gelation in a very different 
way depending on the strength of the binding energy V. 
In particular, in the regime of strong binding close to 
diffusion-limited aggregation, the T dependence of the 
gelation time is governed by Eq. (17) and the gela- 
tion time decreases upon increasing T because the dif- 
fusive transport is enhanced at higher T. In the opposite 
limit of lower V, instead, the gelation time obeys Arrhe- 
nius behavior and increases upon increasing T because 
of the slowing down of aggregation caused by the en- 
hanced thermal breakup of the bonds. Hence, both these 
predicted behaviors appear physically meaningful in the 
two opposite regimes. 

The gelation time as a function of the colloid fraction 
<t> is shown in Fig. 4 for df — 2. Upon increasing the 
attraction, the power-law decay with the exponent —3 — 
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FIG. 4: (color online). Gelation time calculated as a function 
of the colloid fraction cf> for different values of the attraction 
energy V. df = 2. From top to bottom: V — lksT, V = 
2k B T, V = 4fc s T, V = Wk B T. 
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FIG. 5: (color online). Gelation time calculated as a function 
of the attraction energy V for different values of the colloid 
fraction (f>. df — 2. From top to bottom: 0=1- 10~ 4 , 
0=1- 1(T 3 , = 0.01, = 0.1. 



(df/3) predicted in the limit of weak attraction gradually 
decreases and melds into the limiting power-law 0^—1 — 
(df/3) at high attraction. 

The behavior of the gelation time as a function of the 
attraction is plotted in Fig. 5. Three different regimes can 
be identified. At low attraction, V/ksT < 2, the gela- 
tion time decays very rapidly from its asymptotic value 
at V = and melds into a second regime where the de- 
cay with the attraction is exponential. Clearly the first 
regime is dominated by the scaling t g ~ [A/(V / /fcsT)' 3 ] 2 , 
whereas the second regime is controlled by the scaling 
t g ~ [{D/5 2 )e~ v / kBT } 2 , according to Eq.(30). This pre- 
dicted exponential scaling regime is confirmed by ear- 
lier experimental observations 45 . In both these regimes 
the dissociation rate dominates over the association rate, 
i.e. fc_ ^> 4k + ((f)/v). Upon further increasing the at- 



traction, however, the dissociation rate becomes increas- 
ingly smaller in comparison with the association rate as 
4k + (<fi/v) ^$> fc_ and the exponential behavior flattens 
out into a plateau where the gelation time is independent 
of V. This latter regime recovers the diffusion-limited 
irreversible aggregation which is characterized by per- 
manent bonds since dissociation is now infinitely slower 
compared to association. 



IV. STATE DIAGRAM OF ATTRACTIVE 
COLLOIDAL MATTER 

It is possible to summarize the model predictions in a 
state diagram of attractive colloids. In order to be con- 
sistent with previous studies, we introduce the effective 
temperature r customarily defined as^* 



1 
12 



exp 



-V 



(31) 



where a = 2a is the colloid diameter. Clearly, low values 
of r correspond to high attraction energies, and there is 
also a close relation with our dissociation rate, with low 
values of r corresponding to low values of the dissociation 
rate fc_. The state diagram is plotted in Fig. 6. Contin- 
uous lines are calculated for specified values of t g and 
refer to nonequilibrium boundaries. In practice, every 
continuous line separates systems that undergo gelation 
on a time scale t < t g , at higher <fi (i.e. to the right of 
the curve), from systems that undergo gelation at t > t g , 
at lower cf) (i.e. to the left of the curve). All curves are 
plotted for df = 1 since most experimental observations 
of gelation in this regime report values of df close to this 
value, at least for (f> < 0.2. We have checked that chang- 
ing df in the range 1.5 -j- 2.5 does not alter the results 
qualitatively. From these curves it is evident that at least 
in the lower half of the state diagram the gelation process 
is a strongly nonequilibrium process and the transition 
from sol to gel depends crucially on the time-scale of the 
experiment. In particular, if the time of the experiment 
is short compared to the gelation time t g , gelation can- 
not be observed and the system appears liquid-like and at 
most composed of freely diffusing clusters. On the other 
hand, if the time of observation is long compared to the 
gelation time, a transition from a liquid-like material into 
a solid-like one will appear. Hence, in light of our results, 
the observation of so called "equilibrium" clusters in the 
absence of gelation in purely attractive colloids 46 might 
have been due to the time scale of the experiment being 
short compared to the theoretical time scale of gelation 
for those conditions (and in fact the attractions reported 
in RefpS He well in the lower regions of our diagram). 
The situation might be different, however, in the case of 
charged colloids where the electrostatic repulsion plays 
a major role giving rise to further effects^ that are not 
considered in our analysis. The role of the time coor- 
dinate on the gelation transition has been neglected in 




FIG. 6: (color online). State diagram of attractive colloids. 
Continuous lines represent gelation lines (see text for expla- 



nation). From left to right: t g 
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= 0.002, tg 
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0.01, 
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t g — 1 ■ 10 . MR (dashed) line: this line represents the set of 
points on the diagram for which the early-time microscopic 
reversibility condition k+N = fc_ holds exactly. B (dotted) 
line: binodal line. S (dotted) line: spinodal line. 



many previous studies of colloidal gelation, both exper- 
imental and computational, despite being a key control 
parameter in all nonequilibrium transitions. 

In the state diagram we have also plotted the early- 
time microscopic reversibility condition (DB line in Fig. 6) 
represented by the equality between the rates of the mi- 
croscopic forward and backward process, in this case as- 
sociation and dissociation, respectively. The early-time 
microscopic reversibility condition is expressed by the 
equality: k + N = fc_. Regions of the state diagram ly- 
ing much below this line are certainly away from ther- 
modynamic equilibrium and in fact it is below the line 
that we observe the most striking dependence on time. 
Above this line, the gelation lines at different t g tend to 
become more closely spaced together until they almost 
merge together in the top part of the diagram. Inter- 
estingly, the region where all gelation lines (correspond- 
ing to gelation times separated by up to 10 orders of 
magnitude) practically merge, coincides with the binodal 
(liquid-liquid) phase-separation line calculated here using 
a simple mean- field Bragg- Williams approach. Hence, 
from the point of view of kinetic theory, gelation in the 
upper region of the diagram above the binodal line, is an 
extremely unlikely occurrence even over extremely long 
times. From the point of view of equilibrium thermody- 
namics, no gelation can happen above the binodal be- 
cause here the free energy of the system is controlled by 
the unfavorable entropy of mixing which favors the sol 
state (that is obviously more "mixed up" than any ag- 
gregated state) and makes macroscopic aggregation ther- 
modynamically forbidden. 



The metastable region: nucleation and 
liquid-liquid phase separation 



In Fig. 6 we have also plotted the spinodal line and the 
metastable region (delimited by the binodal and by the 
spinodal) deserves careful analysis. Although we have 
plotted our nonequilibrium gelation lines also in this re- 
gion, they should be taken as purely indicative because 
in the metastable region, gelation is replaced by nucle- 
ation leading to liquid-liquid phase separation. Indeed, 
the metastable region appears to be centered upon the 
early-time microscopic reversibility line which is a nec- 
essary condition for the system to be close to the mi- 
croscopic equilibrium between association and dissoci- 
ation and for detailed balance to be satisfied (see also 
ZeldovictPS for the detailed balance principle within the 
context of phase separation). Under these conditions, 
Eq.(l) leads straightforward to homogeneous nucleation, 
as shown in the Appendix. Nucleation leads to the for- 
mation of compact df = 3 clusters which can be seen 
with a simple calculation. 

Being close to the microscopic reversibility line where 
detailed balance can be satisfied (see the Appendix), we 
can write down the free energy for the formation of a 
nucleus or cluster. As in nucleation theory, the nucleus 
is treated as a macroscopic object which allows us to for- 
mulate the free energy of a single cluster. The Gibbs 
free energy contains two contributions. One is the vol- 
ume enthalpy arising from the bonds that are formed: 
q> v = -( Z k/2)V = -(z(R k /a) d f/2)V where k denotes 
the number of particles in the cluster and V the bond 
energy, as usual, and z is the mean number of near- 
est neighbors. The other term is the energy spent to 
create the interface between the cluster and the solvent. 
For a fractal the interface is intrinscally discrete and the 
effective surface can be estimated as the surface occu- 
pied by the particles in the outmost shell of the cluster: 
£Tra 2 dfk( df ~ 1 " df , where we used that the number of par- 
ticles in the outermost shell of a fractal cluster is equal 
to dfk( d} ~ r >' d f . In the limit df — 3 one recovers inR 2 , 
for the cluster surface. 

Therefore, the surface energy is equal to: $ s = 
4n r ya 2 dfk^ df ~ 1 ^ df , where 7 is the surface tension. The 
free energy of the cluster is given by: $ = $„ + 
$ s = -(z(R k /a) d f/2) +A-K^a 2 dfk^ d i- V >l d f. The cluster 
growth at equilibrium occurs along the path of minimal 
free energy. For all values of the parameters involved in 
the free energy, minimization of $ with respect to df in 
3D gives df = 3. We did not consider the <i/-dependcncc 
of z. However, since z is a monotonically increasing func- 
tion of df (because the average density of a cluster in- 
creases upon increasing df at a fixed size) , it is clear that 
accounting for its dependence on df leaves this result un- 
changed. 



B. The unstable region: nonequilibrium 
aggregation 
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V. COMPARISON WITH EXPERIMENTS 



A. Complex viscosity 



Hence we have shown that close to thermodynamic 
equilibrium the clusters are compact objects with df = 3 
which is in agreement with many experimental and sim- 
ulation results presented in the past. As a consequence 
in the metastable region in between the binodal and the 
spinodal line the system is more likely to undergo ag- 
gregation into compact aggregates with the growth law 
R ~ (Dt) 1 / 3 derived in section II. b. 3, and there is no 
competition with gelation. 

It is possible, however, that a solid-like state is formed 
following spinodal decomposition if the volume fraction 
in the dense phase reaches the critical volume fraction 
for the attractive glass transition predicted by Mode- 
Coupling theories, according to a mechan ism tha t has 
been recently discussed in several studies^ 12 * 49 ! Al- 
though this mechanism is certainly a good candidate to 
explain a fluid-solid transition in the proximity of the 
spinodal line, its application is however limited to the 
region of the phase diagram close to the microscopic re- 
versibility line. For deeper quenches well below this line, 
we have £;_ <C k + N and the time scale associated with 
the rearrangement of the bonds on the way to equilibrium 
(i.e. along the minimization of $) is thus longer than 
the time scale on which new particles join the cluster to 
form new bonds. The incoming particles stick irreversibly 
onto the particles in the outer layers of the clusters which 
leads to df < 3 according to the well-known mechanism 
of diffusion-limited fractal aggregation. 

Summarizing, the emerging picture is that sufficiently 
close to the microscopic reversibility line in Fig. 6, the 
clusters can minimize their free energy which leads to 
df = 3. With this input, our analytical solution to the 
master kinetic equation for aggregation gives the growth 
law R ~ t 1 ' 3 , in agreement with the coarsening kinetics 
of spinodally decomposing systems^. Upon departing 
from the equilibrium conditions towards higher attrac- 
tions, the relaxation time over which the cluster mini- 
mizes its free energy becomes long compared to the time 
scale of bond- formation due to the microscopic imbalance 
between association and dissociation. As a consequence, 
df < 3 because the energy minimization cannot be com- 
pleted on the gelation time scale and aggregation pro- 
ceeds with a faster kinetics given by R ~ t x ' dl according 
to Eq.(19). Hence, under these conditions deep inside 
the unstable region in the state diagram, gelation cannot 
be driven by spinodal decomposition because the kinet- 
ics associated with the spinodal-like coarsening is slower 
than the kinetics associated with nonequilibrium fractal 
aggregation. As we are going to see below, this is already 
the case with a relatively mild attraction energy such as 
12k B T. 



The model for the steady shear viscosity of the sys- 
tem as a function of the attraction presented in section 
II. b. 4 can be used within a generalized hydrodynamics 
approach which bridges the gap between the hydrody- 
namic (small u) and the kinetic (large uj) regimes and 
thus provides predictions for the rheological response in 
the whole frequency spectrum. According to generalized 
hydrodynamics, the constitutive relation is written in the 
following Maxwell forrrP^ 



1 d 
G^dt 



d ( dr x 



>xy 



dr. 



dt \ dy dx 



(32) 



where r x and r y are the x and y components, respectively, 
of the microscopic displacement field, and <j xy is the shear 
stress. Upon Laplace-transforming the above equation 
we obtain the complex viscosity as: 



rf{w) = 



Gr, 



— iu + 1/tm 



(33) 



where Goo represents the instantaneous (u —> oo) shear 
modulus and tm = v/Goo is the Maxwell relaxation time. 

Below the transition where the system is fluid, we use 
experimentally measured values of Goo reported in pre- 
vious studies ^ for all calculations. As -q — > oo, Goo is 
the one of the solidified system and can be evaluated as 
the affine contribution to the shear modulus of a disor- 
dered lattice of harmonically-bonded particles. The lat- 
ter, consistent with our picture, are the clusters present 
in the system at the time at which <f> = 0.6422 j n 
fact, whilst the low-w shear modulus is strongly affected 
by nonaffinity^, the high-w one reduces to the affine 
shear modulus given by the following mesoscopic the- 
ory^ Q x = Q A — j^N c z c n<7 2 . All the parameters in 
this formula can be evaluated as follows. The harmonic 
spring constant is k ~ V/S 2 , where V is given by Eq.(35 1, 
and S ~ 10 nm is the hydrophobic attraction range. The 
lattice constant is given by the average cluster diameter. 
This can be calculated using the CMD Eq.(3) evaluated 
at the time at which = 0.64, and the result is a c ~ 2\im. 
N c is the number density of clusters at the random close 
packing and is given by N c — </> c /(7rcr?/6). Finally, the 
coordination number at the random close packing has to 
be^ z c ~ 6. Hence, there are no nontrivial adjustable 
parameters. 

The observable quantity which is measured in the ex- 
periments is the modulus of the complex viscosity, de- 
fined as \rf | = \Jr\i 2 + rjii 2 which then gives: 



W\ = 



Gooh 



Goo/v 2 +u 2 



Gooh 2 + u 2 /G c 



(34) 
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FIG. 7: (color online). Model calculation of the complex vis- 
cosity of thermosensitive colloids using the theory presented 
here [Eq.(34)]. See text for the estimate of Goo- 4> = 0.126, 
ai = 85 ■ 10 -9 nm, uj — Is -1 . The attraction energy is cal- 
culated using Eq. (31) and T c = 33.3°C. As in the experi- 
ment, the temperature is a ramp function of time according 
to T = 30.0+ (i/342.85) °C. From right to left: V = 0M B T, 
V = l.lk B T, V = 2k B T, V = 6k B T, V = 12fc s T. 



From this expression we can extract the most interesting 
limit which gives |tj*| = Goo/w when r\ = oo, i.e. at the 
fluid-solid critical point. 

Eq.(34), together with the microscopic model for the 



cluster evolution as a function of time and attraction, can 
be used to study the rheological response of the system 
upon varying the microscopic parameters. Figf7] displays 
the results of model calculations. The theory predicts 
a sharp rheological fluid-solid transition as a function of 
time upon varying the attraction V(T) provided that the 
final attraction strength V°° is larger than a threshold, 
i.e. V™ — 2fceT for the attraction range used here S — 
lOnm, typical of hydrophobic attraction^] For V°° < 
V^ there is no such sharp transition and gelation is a 
very slow process such that the final solid state may not 
be attained at all on the time scale of observation. The 
curves in Fig. 7 refer to df — 2 and we have checked that 
changing df between two and three has basically no effect 
on the macroscopic transition, apart from shifting the 
solid plateau to somewhat higher values upon increasing 
df, because R c at dynamical arrest clearly is smaller for 
larger df and G m ~ 1/Rq. 



B. Experimental system and interaction 

The colloidal particles used here are the same ones 
characterized in detail in RefpS. They consist of a 52 
nm radius solid polystyrene core onto which a polymeric 
network of crosslinked PNIPAM with T-dependent thick- 
ness (~ 50nm at T < 32°C and ~ 33nm at T > 32°C) 
is affixed. All electrostatic interactions are fully screened 
by the addition of 5.10 -2 M potassium chloride. Attrac- 



FIG. 8: (color online). Comparison between the complex vis- 
cosity from the theoretical predictions of the model (lines) 
and the complex viscosity measured experimentally (sym- 
bols). Volume fraction: (A) 0.126, (□) 0.168. Solid line: 
theoretical prediction using the model presented here with 
= 0.126. 



tion is induced by the hydrophobic effect as the particle 
shell upon collapsing becomes hydrophobic upon increas- 
ing T. The bonds due to hydrophobic attraction be- 
come increasingly long-lived upon increasing T because 
the rate of bond dissociation fe_ decreases with V, ac- 
cording to Eq.(15). As for two-level systems, V is a sig- 
moidal (Fermi-type) function of T and goes from zero up 
to a plateau value V°° across T c according to the follow- 
ing equation derived in RefP^ 



T/oo 



1 



l+cxp[-{T c -T)AS/k B T] 



(35) 



where AS ~ — 1420 A; s is the entropy change across the 
transitiorPS. The latter parameter is an experimental 
input which is fixed by the surface chemistry of the par- 
ticles used in the experiments and it controls the rate of 
variation of V with T. The behavior of Eq.(35) for the 
experimental system used here is shown in Fig.l. 

The complex viscosity as a function of T was measured 
at u = Is -1 in a stress-controlled rotational rheomcter 
MCR 301 (Anton Paar) where T is varied using Peltier 
elements with accuracy of ±0.1°C. The modulus of the 
complex viscosity \rj*\ was measured at a strain 7 = 0.01 
upon varying T at a rate of 0.175°C/min. The solid 
fractions of the dispersions investigated are in the range 
6 — 12% which corresponds to initial volume fractions 
occupied by the particles before the onset of the gelation 
cb<0.26. 



C. Comparison 

In the theoretical calculation we take df — 2 in agree- 
ment with previous studies^. The comparison between 
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the theory and the experiments is shown in FigjHJ It 
is seen that the complex viscosity is constant up to 
T ~ 33°C which is the T value at which the attractive 
interaction sets irPSl. Within this regime, the response, 
according to our viscoelastic model is completely dom- 
inated by the dissipative part, \rj*\ f=a r\. The theory is 
able to capture the very sharp jump of \t)*\ which grows 
by several orders of magnitude within a fraction of degree 
Kelvin shortly after the onset of attraction. Thus, given 
the fact that no nontrivial adjustable parameter has been 
induced, the agreement of theory and experiment can be 
regarded as excellent. 



VI. CONCLUSION 

In conclusion, we proposed a theoretical, fully analyt- 
ical model that bridges the microscopic physics of col- 
loidal interactions and physical bond-formation and dis- 
sociation, with the macroscopic rheology and the col- 
loidal gelation transition. According to the model pre- 
dictions, for shallow attraction energy V, growth oc- 
curs with a lag time 9 ~ De~ v ^ ksT . The size grows 
as R ~ [1 + (t/9)] 1 ^ df , asymptotically recovering the 
Lifshitz-Slyozov law R ~ (Dt) 1 ' 3 in the equilibrium limit 
df = 3. The gelation time t g decays as a power law 
of the volume fraction with an exponent which is both 
attraction-dependent and ^-dependent. In the limit of 
weak attraction, t g ~ 0- 3 ~( d // 3 ), whilst t g ~ ^-M^z/a) 
is found in the infinite attraction limit. The theory 
also yields a state diagram of attractive colloids. In 
the metastable region, an equilibrium-like nucleation sce- 
nario can be recovered as long as detailed balance holds, 
at least approximately. Although it is possible that gela- 
tion is driven by spinodal decomposition close to the 
spinodal line, identifying a gelation boundary which is 
independent of time for deeper quenches seems unfea- 
sible and nonequilibrium aggregation with df < 3 is a 
faster process than spinodal coarsening with df — 3. Our 
results are supported by experiments on a model suspen- 
sion of thermosensitive colloids in which V is increased 
from zero in the same sample by raising T. The compari- 
son shows that gelation already for V — 12ksT is a sharp 
nonequilibrium fluid-solid transition with df = 2, and 
culminates with a solid random packing of clusters with 
infinite viscosity. The solidification is irreversible unless 
one switches off the attraction by reverting T, which is a 
unique feature of the experimental system under study. 

We have also presented experimental observations of 
colloidal gelation in a system of thermosensitive colloidal 
particles where the attraction is varied from zero up to a 
maximum in the same system by simply varying T. The 
complex viscosity as a function of T can be quantitatively 
modeled with the theory presented here and exhibits a 
sharp gelation transition with df = 2. The comparison 
between theory and experiments also indicates that al- 
ready with a mild attraction energy of 12kgT short-range 
attractive colloids undergo fast nonequilibrium gelation 



with df — 2 and with no detectable hallmarks of spin- 
odal decomposition. Our analysis suggests that the latter 
mechanism is more likely to play a role and affect gelation 
for weaker attractions close to the spinodal line where the 
rates of two-particle association and dissociation are close 
in value and aggregation occurs with df — 3. 

Clearly, the framework presented here opens the 
unprecedented possibility of engineering functional 
materials that turn into solid on a desired time-scale and 
that can be switched in a fully controlled and reversible 
way between fluid and solid states. 

APPENDIX: MICROSCOPIC REVERSIBIL- 
ITY, DETAILED BALANCE, AND THE ON- 
SET OF NONEQUILIBRIUM AGGREGATION 

Let us consider one-step association and dissociation pro- 
cesses by which one particle joins a cluster and disso- 
ciates from a cluster, respectively. According to the 
definition given by TolmarPD the principle of micro- 
scopic reversibility is satisfied when for every cluster 
j the absolute frequencies of particle attachment (to a 
j — 1 cluster) is equal to the absolute frequency of par- 



ticle detachment (from a j cluster): kK 



k J , here 



n represents the number of monomers in the system. 
Clearly, in the early stage of the aggregation process, 
if the latter is sufficiently slow, n can be taken as a con- 
stant and equal to the initial concentration of monomers 
N . Then, upon neglecting cluster size dependencies, we 
recover what we called the early-time microscopic re- 
versibility condition k + N = fc_. When this condition 
is strongly violated and the association process to form 
dimers is initially much faster than the dimer dissocia- 
tion process, k+N ^> £;_, then the formation of trimers 
etc. becomes a faster process. As a consequence of the 
fast trimer formation, the equilibrium concentration of 
dimers given by n^ 9 = (N — n(t — > oo))/2 where n(t) is 
given by Eq. (7), cannot be reached. Then, the condi- 
tion k+(n eq ) 2 — fc_(n2 <3 ) 2 for j = 2 cannot be satisfied. 
In turn, this fact implies that also the detailed balance 
condition kP + ~ {n^) 2 = k^(n^ q ) 2 which has to be satis- 
fied V j, cannot hold since it is already violated at least 
for j = 2. In view of this argument, it is clear that the 
condition k + N = fc_ that we encountered multiple times 
and plays a crucial role in the development of our theory, 
is tightly connected with the detailed balance principle 
and it may be regarded as a necessary (although prob- 
ably not sufficient) condition for detailed balance to be 
satisfied. Hence, this argument also provides a simple 
but robust quantitative criterion to assess whether the 
aggregation process is nonequilibrium or not, with the 
inequality k + N > fc_ marking the crossover from the 
equilibrium aggregation (nucleation) dynamics into the 
nonequilibrium one. 

Finally, we should interpret the approximation under- 
lying Eq.(9) within this framework. The truncation to 
first-order in time is equivalent to neglecting the kinetic 
equilibrium plateau for the dimers and to restricting the 
analysis to the initial growth stage only which occurs 
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before the plateau. The time at which the plateau is 
reached can be easily estimated from Eq.(7). Since the 
tanh reaches a plateau when its argument is ~ 2, this 
gives t eq = Ayjk-{k- + Ak+N) for the time required to 
establish the law-of-mass action kinetic equilibrium for 
the dimcrs. If the arrival frequency of the third particles 
is k + N ', when k + N > t eq then only the linearly grow- 
ing part of a(t) from Eqs.(7)-(8) is necessary to arrive at 
K c ff . Then the inequality to be satisfied for our approxi- 
mation to be fully justified is k + N > 4-^/fcI (fc_ + Ak + N) 
which can be rewritten as k + N > 16fc_. Hence, any re- 
sult presented here in the regime k + N < 16k- N should 
be regarded as merely illustrative. This fact is however 
not worrisome because in that regime the crossover into 
nucleation takes place and our kinetic theory anyways 
has to be replaced by nucleation and coarsening at some 
point. 

APPENDIX: NUCLEATION KINETICS 
In this Appendix we shall see how the classical homoge- 
neous nucleation is recovered when the detailed balance 
condition is satisfied and the aggregation occurs with 
df = 3 by minimizing the free energy of clustering <I> 
according to the mechanism discussed in section IV. a. 
Recall that the master kinetic equation for aggregation 
and breakup reads as: 



dN k 1 



i-\-j — k 



V K+.NiN, - N k V K+Ni 



K-N k + £ K7 k N t 

i=fe+l 



(36) 



Let us assume that detailed balance is satisfied such that 
aggregation occurs under equilibrium conditions by min- 
imizing the free energy $ = $ v + $ s = -(z(R k /a) d f /2) + 
A-K r )a 2 dfk^ df ^ 1 >' d! . Minimization gives df — 3, mainly 
to maximize the mean coordination number z. Then only 
detachments of individual colloidal particles are relevant 
whereas body-fragmentation plays no role because it in- 
volves the breakup of many bonds, which is a much slower 
process. Furthermore, clusters form as a result of spon- 
taneous fluctuations and the growth of each cluster is 
independent of the behavior of the others. Clearly, un- 
der these conditions, the above equation can be rewritten 
accounting only for one-step association and dissociation 
processes: 

r ' 4 '" K+iN k -K-N k ^j^ k+11 iM k+1 . 

(37) 

where we have incorporated the constant factor N\ = N 
inside the association rate constants, as the monomer 
concentration is assumed to be constant for the slow pro- 
cess. To shorten the notation we put K k _ 1 x = Kt_ 1 etc. 
and rewrite this equation as: 

^ = i^ATfc-i - K k N k - K-N k + K- +l N k+1 . (38) 

Since the attraction is relatively weak and thermal dis- 
sociation is important, the principle of detailed balance 



dt 



*T-i,i^-i 



is applicable in this limit. Hence, we now introduce the 
equilibrium or steady-state concentration of aggregates 
of size k N k q which is a Boltzmann function of the min- 
imum available work (free energy) $ needed to form an 
aggregate of size k: N^ q ~ exp(— $/fcT). Upon applying 
the principle of detailed balance we have: 



A^i^-i = K q K k 



(39) 



These relations allow us to eliminate from Eq. (38) the 
quantities K k and K k+1 and Eq. (38) becomes: 



dN k 
dt 



K k 



^iJVT/^J 



■K, 



fe-i 



■N k + N k 



(40) 



= K k Nl«[N k+1 /NH 1 -N k /Nl 



(.:<;] 



- K k ^Nl\ [N k /N 



'f 



N k 



-i/K q - 



Let us now transform the discrete distribution N k into a 
continuous one N(x) where x ~ R is a continuous vari- 
able expressing the cluster size. Denoting by A the spac- 
ing along the x axis between the neighboring sizes k and 
k + 1, we have N k = XN(x), N k+ i = XN(x + A) etc. 
Then we get: 



dNk 
dt 



. K{x)N eq {x) [N(x + X)/N eq {x + A) - N(x)/N eq {x)} 

- K(x - X)N eq (x - X)[N(x)/N eq (x) 

-N(x-X)/N eq (x-X)}. 

(41) 
Since A is constant and N, N eq and K vary little within 
the length A, one can do an expansion in power series of A 
up to the first non- vanishing terms, under the assumption 
that N varies slowly and hence N > XN' > A 2 TV". The 
Taylor expansion is done on the terms in (x + A) and 
(x — A) and is centered on x: 



dN(x) 
dt 



K(x)N eq (x) 



K{x-X)N eq {x-X) 



d N(x) ' 
dxN e i(x)_ 

d N(x-X) 



A 



dxN e< i{x- X) 
We now expand the second term on the r.h.s.: 

d N(x-X) 



(42) 



K(x-X)N eq (x-X) 
X 



K{x)N eq {x) 
d 



~dxN ec i(x- A) 
d N(x] 



X-Z- { K(x)N eq (x 
ox 



dxN e i{x) 

d N(x) 



(43) 



A- 



dx N e i(x) 



Upon replacing Eq.(43) in Eq. (42), cancelation of terms 
leads to: 



^-»'!w^<.. 



d N(x) 



dx N e i(x) 



(44) 
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The product in the braces can be rewritten as: 

kn^ r d ( N ^ 

.dN(x 



dx \N e i(x) 



= K- 



dx 



KN(x 



d\nN ec i(x) 
dx 



(45) 



The equilibrium distribution is given by the Boltzmann 
form: N eq (x) ~ e -*(z)/fcT\ Inserting this and upon fi- 
nally replacing in Eq.(44) we obtain: 



dN(x) 



iim + ^-mi <*. 



dx \ dx 



Bt 

where D c g = X 2 K is an effective diffusion coefficient in 



size space. Thus, we obtained a diffusion equation in size 
space for the growth kinetics. This equation can solved 
with Kramers saddle-point approximation to recover the 
standard nucleation rate for systems close to equilibrium: 






$(x*) 
' k B T 



(47) 



where x* denotes the critical nucleus corresponding to 
the maximum of <&, i.e. the size at which the surface 
term and the enthalpy term in $ balance. 

* az302@cam.ac.uk 



1 H.J. Herrmann, Physics Reports 136,153 (1986). 

2 W.C.K. Poon and M.D. Haw, Adv. Colloid Interface Sci. 
73, 71 (1997). 

3 V. Trappe and P. Sandkuehler, Curr. Opin. Colloid Inter- 
face Sci. 8, 494 (2004). 

4 S. Romer, F. Scheffold and P. Schurtenberger, Phys. Rev. 
Lett. 85, 4890 (2000). 

5 P. N. Segre, V. Prasad, A. B. Schofield, and D. A. Weitz, 
Phys. Rev. Lett. 86, 6042 (2001). 

6 V. Trappe et al., Nature 411, 772 (2001). 

7 E. Zaccarelli, J. Phys.: Condensed Matter 19, 323101 
(2007). 

8 P. J. Lu et al. Nature 453, 499 (2008). 

9 S. Buzzaccaro, R. Rusconi, and R. Piazza, Phys. Rev. Lett. 
99, 098301 (2008). 

10 M. J. Solomon and P. Varadan, Phys. Rev. E 63, 051402. 

11 H. Tanaka, S. Jabbari-Farouji, J. Meunier, and D. Bonn, 
Phys. Rev. E 71, 021402 (2005). 

12 T. Gibaud and P. Schurtenberger, J. Phys.: Condensed 
Matter 21, 322201 (2009). 

13 A.P.R. Eberle, N.J. Wagner, R. Castaneda-Priego, Phys. 
Rev. Lett. 106, 105704 (2011). 

14 W. Y. Shih, I. A. Aksay, and R. Kikuchi, Phys. Rev. A 36, 
5015 (1987). 

15 M.D. Haw, M. Sievwright, W.C.K. Poon, P.N. Pusey, Adv. 
Colloid Interface Sci. 61, 1 (1995). 

16 M.A. Miller and D. Frenkel, Phys. Rev. Lett. 90, 135702 
(2003). 

17 S.A. Shah, Y.L. Chen, K.S. Schweizer, and C.F. Zukoski, 
J. Chem. Phys. 119, 8747 (2003). 

18 J. Bergenholtz, W. C. K. Poon and M. Fuchs, Langmuir 
19, 4493 (2003). 

19 G. Foffi, C. De Michele, F. Sciortino, and P. Tartaglia, 
Phys. Rev. Lett. 94, 078301 (2005). 

20 E. Del Gado and W. Kob, Phys. Rev. Lett. 98, 028303 
(2007). 

21 K. Kroy, M. E. Cates, and W. C. K. Poon, Phys. Rev. 
Lett. 92, 148302 (2004). 

22 H.H. Winter and F. Chambon, J. Rheol. 30, 367 (1986). 

23 E. Del Gado, A. Fierro, L. de Arcangelis, and A. Coniglio, 
Europhys. Lett. 63, 1 (2003). 

24 A. Zaccone, H. Wu and E. Del Gado, Phys. Rev. Lett. 
103, 208301 (2009); M. Laurati et al., J. Chem. Phys. 



130, 134907 (2009). 

25 A. Zaccone, J. J. Crassous, B. Beri and M. Ballauff, Phys. 
Rev. Lett. 107, 168303 (2011). 

26 V. I. Agoshkov, P. B. Dubovski, V. P. Shutyaev, Methods 
for Solving Mathematical Physics Problems (Cambridge 
Int. Science Publishing, Cambridge, 2006). 

27 H. Holthoff et al, Langmuir 12, 5541 (1996). 

28 A. Zaccone and E. M. Terentjev, Phys. Rev. Lett. 108, 
038302 (2012). 

29 C. Schneider, A. Jusufi, R. Farina, P. Pincus, M. Tirrell, 
and M. Ballauff, Phys. Rev. E 82, 011401 (2010). 

30 H. A. Kramers, Physica 7, 284 (1940). 

31 R.P. Tolman, The Principles of Statistical Mechanics, (Ox- 
ford University Press, Oxford, 1938). 

32 M. Abramowitz and I. A. Stegun, Handbook of Mathemat- 
ical Functions with Formulas, Graphs, and Mathematical 
Tables (Dover, New York, 1972). 

33 D. Asnaghi, M. Carpineti, M. Giglio, and M. Sozzi, Phys. 
Rev. A 45, 1018 (1992). 

34 E.M. Lifshitz and L.P. Pitaevskii, Physical Kinetics (Perg- 
amon Press, OPxford, 1981). 

35 P. Sandkuehler, J. Sefcik, and M. Morbidelli, J. Phys. 
Chem. B 108, 20105-20121 (2004). 

36 A. Zaccone, D. Gentili, H. Wu, M. Morbidelli, E. Del Gado, 
Phys. Rev. Lett. 106, 138301 (2011) 

37 P. G. Debenedetti, Metastable Liquids (Princeton Univer- 
sity Press, Princeton NJ, 1996). 

38 P. G. de Gennes, Scaling Concepts in Polymer Physics 
(Cornell Univ. Press, Ithaca NY, 1979) p. 176. 

39 L.D. Landau and E.M. Lifshitz, Fluid Mechanics 
(Butterworth-Heinemann, Oxford, 1987). 

40 C.I. Mendoza and I. Santamara-Holek, J. Chem. Phys. 130, 
044904 (2009). 

41 W. Schaertl and H. Sillescu, J. Stat. Phys. 77, 1007, (1994). 

42 G. Lois, J. Blawzdziewicz, and C. S. O'Hern, Phys. Rev. 
Lett. 100, 028001 (2008). 

43 W. Goetze and L. Sjoegren, Z. Phys. B 65, 415 (1987). 

44 Z. Cheng, J. Zhu, P. M. Chaikin, S.-E. Phan, and W. B. 
Russel, Phys. Rev. E 65, 041405 (2002). 

45 W.C.K. Poon , et al. Faraday Discuss. 112, 143 (1999). 

46 P.J. Lu et al., Phys. Rev. Lett. 96, 028306 (2006). 

47 A. Stradner et al., Nature 432, 492 (2004). 

48 Ya.B. Zeldovich, J. Phys. USSR 12, 525 (1942). 



15 

49 S. Manley, et al. Phys. Rev. Lett. 95, 238302 (2005). 53 C. S. O'Hern, et al. Phys. Rev. E 68, 011306 (2003). 

50 J. P. Hansen and I. R. McDonald, Theory of Simple Liquids 54 A. Fernandez-Nieves et al, Langmuir 17, 1841 (2001); M. 
(Academic Press, London, 2006). Rasmusson, A. Routh, and B. Vincent, Langmuir 20, 3536 

51 I. Deike, M. Ballauff, N. Willenbacher, and A. Weiss, J. (2004); M. Ballauff and Y. Lu, Progr. Polym. Sci. 36, 767 
Rheol. 45, 709 (2001) (2011). 

52 A. Zaccone and E. Scossa- Romano, Phys. Rev. B 83, 
184205 (2011). 



